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det The purpose ose of this ps paper is to asain the: nature ire of unconfined flo flow to single 
and multiple we Ww ells, and to present a ‘method of solving pre oblems associated w vith 
3 this typ pe of flow. ‘The effect of the ¢ capillary fringe on the location | of the fr ee 
= and the form of the flow patterns, the zone of validity of the Dupuit 
s equation, the shape | of the fr ee ‘surface near the well, and the variation in the 
"stream s surface ‘spacing ar are all discussed. A functional | relationship independ- E. 


of radius | of influence i is hed, vials the variables at the Ww well; 


this relationship applies to both single and wells. fund..me di- 


-mensionless | parameter consisting of a ratio of Froude’s Rey eynolds” numbers 


is formulated 1 that characterizes the : shape of the cone of depression on around a a 


Ww ell. The concepts of w ell efficiency and effectiveness are clarified and guides 


are d for their correct use. 


presente elr correct use. 
% 
_ INTRODUCTION 
y Notation. —The letter symbols introduced i in this paper are 2 defined w here 


they first appear ‘and are assembled alphabetically i in the Appendix : for con- “2 


7" The need for a clearer understanding « of the flow w to unconfined w wells and for 
better er methods a of solving. the problems ari ising , from this flow has arisen asa 


‘result. of extensive pumping of ground w ater and the demand for a more eco-— 


-nomical use of the w ater. The investigation reported in this paper was under-| 
‘taken in an effort. to contribute further knowledge both the nature of 


flow and the method of solution. 


gs Flow into’ wells ¢ ean ‘be divided into tw o broad classifications, depending 
the boundary conditions. The first classification is confined flow, 
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—lTrrig. Engr. Trrig. Div., SCS; and Asst. Prof. of en Utah State Agri. College, Logan, Utah. oa 


Proceedings—Vol.. 78 Au ust, 1952 | o. 142 — 
si 
BY Vay ‘MASCE 
| 
= 
| 
— 
| | 
b 
— 
— 
— 
| 
— 
— 


M MULTIPLE WE: LLS 


‘Since confined flow occurs between fixed is susceptible 
rather complete analysis based on re relatively few approximations and assump-— 
: tions, with the result that the theoretical solutions conform very closely to the 
observed conditions. Consequently, considerable knowledge is available re- _ 
garding the theory and practical aspects of such flow. 
 Unconfined flow is much more involved, however, because of the existence 
of an unbounded water surface. _ Since the shape of the surface i is not known, 
an ay approximate . analysis based on extensive assumptions m must be made, — 
> _ imposing rather severe limitations on the resulting solution. — ¥ The general ace 
7 ceptance of the solution, | without understanding the inherent limitations, has 7 
=. produced a a certain n degree of complacency s among many | practical ground-w -water i 
"engineers, as w well as sa tendency t to disregard any observations at variance with 7 
the expected solution. - It is important, also, that the theoretician clearly ap-— 
_ Praise the limitations - imposed by nature on the too extensive application of 7 


4 equations of ground- water flow; nature rarely provides the idealized case as- 


sumed in most theoretical developments. However, the s simplified situation | 


. must remain the basis of study, with the knowledge that it does not occur, and : 


the results must be modified accordingly. 
In this study, the usual assumptions of steady flow ‘ through a homogeneous ~ 


7 _ porous medium hav have | been made. _ For the confined flow, the bounding surfaces 


a are assumed to be parallel, horizontal, imper ‘meable layers; for the unconfined 
flow, the wells are assumed to penetrate the medium 7 


: overlying an impermeable, horizontal 


A analysis of flow through porous has been | extensively 


covered d by several writers.” "3 The essentials of these analyses, as applied to 
- the case of stes ady flow into wells, will be summarized in order to interpret the 
Significance of the experimental results. 


hich v is velocity, K is the of h is the 
head, and s is the length along a stream line. Eq. 1 is combined with 
partial differential equation of continuity | Vv=0 for confined flow, the La- 
place equation (V? h =| 0) is” obtained. W hen 1 the boundary conditions for 


— confined radial flow are applied and the differential equation is integrated, the 
following eque etion for the shape of the aan surface ofa 0 artesian ; 


_2“Flow of Ground Water,” by Cc. E. Jacob, “Enginecring John ‘Wiley Sons, Ine., » New 
ork, N. Y., 1950, Chapter 5. 


8“The Flow of Fluids Through Media, “i M. Muskat, J. LW. Baw Edwards, Ann 
Arbor, Mich., 1946. 


a which the water is restric her impermea 
= 
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— 
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in which Q is the total discharge, t¢ is the ihn, and r is the radius of the | 
= . Since the Laplace ce equation is I linear i in h, the solution for multiple = 


is simply the sum of the solutions for each individual well. _ 
sy As was previously stated, the existence of a free surface greatly — ow 


the analysis. One writer has gone so far as to claim that the problem cannot 
be solved with the present knowledge of mathematics. — This, of cou course, is an 

_ extreme statement of the complexity of the problem, for through the use of 
simplifying assumptions, s, clearly understood, a solution may be obtained that 
will be of considerable value. The failure to understand these assumptions © 

_ and their e effect upon the solution has been the r root of a great deal of confusion — 

in the minds of many practical a as well as theoretical investigators. wes . 
‘The classical analysis of this problem was presented by Jules _Dupuit in 

, 1863. He assumed horizontal flow throughout a homogeneous s material under- 


well completely penetrating the perme- 


Tain by by an eyo sie stratum with a w 


having a height equal to the ‘depth of the water in well. The 


ofa an analysis based on these assumptions is the formula commonly known as 
the: equation: ran 


K 


24 Limitations « of the Dupuit Equation -—Since the free surface approaches the 

_ horizontal at a considerable distance from the well, it can be seen that the as 

sumption of horizontal flow through vertical potential surfaces Ww ill become m more 
accurate as the radius i increases, , with a , resulting i increase in the : accuracy of nm 
equation. 

_ The surface given by the Dupuit equation intersects the well at the level of 4 

the water in the well, whereas it is a known fact that the free surface intersects — 


Mid 


the well above this point, giving rise t¢ to what i is commonly referred to as a Sapal 


- of seepage that i increases in extent as the drawdown increases. — ~ Moreover as 
the fluid approaches the w ell, the free-surface slope i increases , causing the po- 
- tential si surfaces (that must be normal to the flow surfaces) to depart more and 
from the cylindrical form assumed in the derivation. these reasons 
oe actual free surface would be above that calculated by the Dupuit equation, q 
the difference being ig greatest at the edge of the well. 
Effect of a Capillary Zone—In addition to the flow condition ns already 
mentioned, the ef effect of the e capillary zone m must be: considered. ‘This i is 
true for model studies w herein the capillary 1 rise. be an 
percentage of t the depth of flow. — Ordinarily it is sufficiently accurate to ‘say . 
that the free surface i is the surface below» which the voids are saturated with | - 
fluid. > However, when a capillary zone becomes important because of its 
relative s size, this d definition i is no longer su: sufficient and the free surface can best 
_ be defined as the surface of atmospheri iG | pressure. e. R.D. Wyckoff, H. G. Botset, 


and Mr. Muskat,‘ and others have found that the > potential distribution —— { 


eae. “Flow of Liquids Through Porous Media "Under the Action of Gravity,’’ by R. D. Wyckoff, H. G. 
and M. Muskat, Physics, Vol.3,1932pp.90-114.0 
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across the. atmospheric- -pressure surfs ace into the cap capillary zone. “This potential 
- distribution i in the capillary zone will g giv e rise to additional flow that cannot be 
Several ators, recognizing g the limitations of the Dupuit. equation, 
have conducted m model tests to ascertain closely the true nature of the complex 
problem | of unconfined flow. In 1932, Messrs. WwW yckoff, Botset, and Muskat! 
constructed a 15° sand sector with glass sides and piezometers s connected to i 
pevrnee . They found that the base piezometr ic heads, rather than the free sur-_ 


face, were given by the Dupuit equation. _ By injecting dye into the sand at 
the inflow face and tracing its movement to the well, they found that the capil- 
zone contributed a sizeable discharge >. These measurements were used to 


MAES? 
substantiate a modification of the Dupuit equation that allows for ‘the increased 7 


flow w resulting from capillar ity. Apparently no no attempt wa was made to deter ermine 


the effect of the capillary flow on the free surface. 


E. Babbitt, M. ASCE, and David H. A.M. ASCE, 


undertaken a a very extensive three- phase » study of un unconfined flow. ay The fret 


phase o of the study was a series of electrical-analogy tests on a thin carbon w edge 


representing a single w w well. In addition to verifying the previously mentioned 


limitations of the Dupuit “equation, they devi eloped | an equation for the nin 


The second phase of the i inv estigation was the construction of a 15° cceual 
sector 100 in. 1. long and 16 in. high having a single well located at the ap apex. < 


” ‘Along one of the sides the piezometric head was measured and the free surface — 


determined. _ The discharges Ww ere compared with those given by the > equation 


developed by the electrical- analogy ‘method and w ere found to be in good 


Astudy of multiple v w ells comprised the third phase. 13- ft aquare con- 
_ tainer wa as filled to a depth of 13 in. with ordinary building sand. Ww ater Ww as” 
admitted through screens located at the corners and flowed toward various com- 


_ binations of 5 wells placed hear the center. The discharges were compared — 
Bess those computed from a multiple-w ell analys sis based on certain stain ing 


assumptions. air a nent was s ed. 
sumption fair creement W ecured. 


a. Since the study by Messrs. Babbitt and Caldw ell represents the ener ex 


-haustive investigation of unconfined flow known to the writer, ‘it 
as the starting point for the experiments. reported herein. n. For that reason an 
appraisal 0 of the results of the Babbitt and Caldw ell tests is needed at 1 this aad 
understand clearly the « objectives of this investigation. 
= First of all, with regard to the electrical- analogy muds, it is felt oe the 
resulting free surface i is a good representation of the true. fr ree surface. How- 


ever, ¢ caution should be used near the well, s since the w edge | w ras very thin and 
s or lack of homogeneity as well as improper r electrode contact: 


any small cracks or 
would be more notice eable in this Tegion. . Moreover, the range of variables — 


_ &“The Free Surface Around, and Interference Between, Gravity Wells,’’ by Harold E. Babbitt = = 
Dav id H. Caldwell, Bulletin Series No. 374, Experiment Station, Urbana, IIl., January 7, 1948. 
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studied was rather small. For maximum ex 


and Caldw ell was : presented i in 


1K 
_ in which the subscript e ‘refers to the external boundary, | and the coefficient C, 
is defined by a curve relating the draw down to the radius. — Of extreme interest — 
is the fact that when. + is plotted against : the logar ithm of the quantity (r/r.), 
a straight line results to : a VE alue | of r/ re = 0.05. _ Since the minimum value of 
tft in the electrical- -analogy experini riments was s only 0. 04 4 and the maximum 0. um 0.087, 
the logarithmic relation holds essentially true to the ‘edge « of the well, pial 


i When the approximation (C, = 0. 3 log r/re) is substituted in Eq. 4, the 


following formula results: = 
los 
_ Notew orthy is the fact that the elevation of the free surface i is expreased as a 


_ linear function of the radius 
Whereas the Dupuit equation is not a linear 
Mopet Stupres or UNCONFINED 


Equipment for “Model Studies. .—The sand- -model studies described herein 
were: undertaken at the laboratory of the Iowa Institute of Hydraulic Research, 


ncipally for the purpose of gaining further information about the free atten 


principa 
= the well. A 90° sector containing sand (Fig. 1) was constructed, in which > 


the inflow screen was 120 in. from a central well having a radius of 1.23 in. i. The 


depth of inflow was maintained at approximately 24 in. in 30 in. of sand. 7 Pwo 
_ other wells, each of 1.2-in. radius, were located 36 in. from the apex . Thew fells 
were placed along planes « of symmetry s so that by the p principle of images the full 
"pattern of wells would be in the center of a circular confining boundary. — One- 
of the p planes of symmetry was made of plexiglas, containing 174 piezometer 
potent omg inflow face was formed « of two layers of copper screen sup- 
ported by i x4 4-in. . timber Tibs. — ~The base of the model was constructed of 
concrete and the 1e sides, ¢ exclusive. of the plexiglas, were of } i-in. plyw wood. Two 
layers of copper screen, soldered to 3-in. hardw are cloth, formed the well, face. “ 


WwW ater entered the sand through the screen at the outer periphery and flowe 


‘toward the wells, at which point the water surface w as maintained at the re- 
quired. level by an overflow control. The discharge was measured by w weighing 


water flowing from the w wells. 
th most 


of f the by w washing. To reduce the amount of air, | Ww hich | 
has a considerable and variable effect upon the permeability,® the snail es 

placed under water. Each layer was approximately 1-in. thick and raked 
58, 


6° of Air ir upon the e Permeability of Soils, by y J. E, Christiansen, Soil Science, Vol, 


The formula propose 
pte 
— 
? 
fj 
— 
— 
| ; 
pier 
— 
— 


q 


. 
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_ Horoaghily the horizontal stratification as well as the orientation of the grains. a 

Procedure for Model Studies ——The sequence of operations in the series of 
nthe 

tests was determined by the ¢ criterion that the water level i in the » sand should 


: reduced. This effect could very well be serious near the well at which pc point f 


_ The e water s surface i in the. central we ell w ‘as low ered to 18. 31 in. -and d the reservoir q 
i wi el to approximately. 24 i in., which was maintained during the remainder of the 
tests. _ Since the sand was 30 ‘in. deep, this differential allowed for the expected | 
maximum capillary 1 rise of 6 in. to be completed r near the outer periphery. - The 
_ piezometric heads, together with the discharge, were used in ee when 
an equilibrium condition Was sapproached. 
strength red food coloring w was placed i in all piezometer This. 
dye served two very | useful purposes in addition to functioning as the piezom- — 
4 eter fluid. | The dye. could be made to enter the flow and trace out the stream _ 


lines (Fig. 2) by either placing additional dye in n the tubes or a 


air ir would be the pores. and the thus be 


> 
i 
— 
— 


tubes were given an initial compression, the velocity of the resulting ¢ dye front 
could be determined. The piezometer readings \ were corrected for the increase 
in ‘specific gravity resulting from the coloring. 
Experimental Tests on Model.—Seven complete. tests were made. Th 


a da well, referred to hereafter as well 1 (Figs. 1 and 3), was operated | at four 


1 


—Dye Fronts Our 
levels: 18.3 in., 12.3 in., 6.3 in., and ‘complete drawdown. (All were 
measured above the base of the model.) A typical result is shown in Fig. 
The fifth test set up consisted of 3 wells in a straight line with the central 


: well maintained at full drawdow n and the outside wells at 12.9 in. elevation. 


— 
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The potential | distribution f for this case was obtained in two steps. ea In the first 
; step w ell 1 was at full drawdown and well 3 at 12.9 i in. - ‘The flow pattern thus 
3 determined along the plexiglas boundary was then at right angles to the line of 
Ww ells. The second step was to close well 3 completely and lower the water sur- 


oF Wa ELL Layour For Movex Tests 

face it in Ww well 2 2 to 12, 9i in. Int this Manner the flow y pattern along t the line of wells — 


was determined. It will be observed that it was necessary to raise the w: vater 
sur ‘face over a portion of the flow area, but the tests were ar ranged so that this 


Z fluctuation alwe vays occurred at w vell 3, the e farthest from the line of piezometers = 


In this way Ly it was hoped to minimize the effect of trapped air. 
a ‘s ~The sixth test w ras made w ith 4 w ells i in a square and the fifth well in the | 


center. | Thep patter n consisted of 3 of the wells in a line along a diagonal operat- — 


© = Actual Free Surface 
Capillary Fringe 


in Reservoir Adjusted Free Surface~_ 


Direction 
Lines 
B *Potential Lines’ Water Surface’ 


My “4 Bey at full draw vdow n, and the remaining 2 


was obtained with wells 1 and 3 at full euaiain an w vell 2 a 12.9 in. “The 


_ flow pattern thus obtained w yas at right angles to the line of 3 w ells at full draw- - 

a down .W hen well 2 was lowe yered to full draw vdown and well 3 was raised to 
12. 9 in., the flow pattern along the line of the 3 wells wa s obtained. _ a5 

Td r The last test, number seven, was W) with all all 3 wells at f full draw ‘down, resulting 
1a patter n of 5 wells under full oper ration ; ‘the results ¢ of this test are shown in 


= & = 
| 
— 
ae 
dge of Well. [Well 
— 


“MULTIPLE WELLS 


Analysis of Test | Results. —The discussion of the sand- model tests will be 


—the Dupuit analysis, the capillary zo zone, the free 
surface, and d the boundary conditions at the 


The Dupuit Analysis.— tion has been shown previously to — 

agree with the eeeenines heads measured along the bottom of the permeable — 


/in Reservoir Edge of Well | Well 


Measured Free Surface 


Potential Lines Water — 
—F.Low FOR Test NUMBER 7 
stratum. Additional experimental verification is shown in Fig. 6, in w vhich. the 
: -Dupuit equation i is the best fit curve | for the experimental data obtained from 
tests one to four. WwW ith the experimental tests indicating that the Dupuit - 
“equation accurately g gives the free surface at considerable distances from the __ 
: well and also the piezometr ic head along the entire base, a re- examination of the , 
-Dupuit analysis will be of interest in ascertaining the reason for ‘agreement in 
regions. The basic assumption of horizontal flow through vertical 


tential surfaces _ is. is essentially met at considerable distances from the well. 
Furthern more, since ce the p potential surfaces are normal to the horizontal i impervi- a 
ous bed, the : assumption of horizontal flow through \ vertical potential surfaces - 


is also reasonable in this ‘region. — The flow between the base and an adjacent. - 

flow surface is 

an =2nrrA 

in which q is the discharge iia een flow surfaces. i’ ‘When it is assumed that 
=Ch and | therefore q = 7 CQ, the expression may b be e integrated t to obtain 

Dupuit equation. Since the resulting equation is verified experimentally, 

_ the assumption is evidently sound, and can therefore be used in further defining _ 

‘the pattern of flow. In brief, this analysis means that the normal distance be- ‘ 


tween stream surfaces immediately adjacent to the i imper meable base remains a 7 


constant fraction of the the impervious base the 


Effect « of the Capillary Zone.—The capillary zone had considerable 


_ upon the results, as its relative magnitude would indicate—the inflow depth 
-_= approximately 24 in. . and the average capillary rise being 5 in. sa From | 
Figs, 4 and 5, the potentials can be seen to extend across the free surface ont 7 


— 
— 
4 — 
| 
a 
4 
rs 
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a Ww well into the capillary zone. This, 0 of course, should be expected s since the flow a 
- of the fluid depends on the e gradient and not the absolute magnitude of the 


~ 
the diflerence in 1 permeability that the fringe any 
_ graded granular material, a fringe rather than a definite clear-cut line will exist — 
at the top | of the capillary zone. 
int ‘The Free Surface. —The inflow surface that must exist in all model studies 7 
ad will be considered first. When the fluid enters the permeable r material, hana 
‘capillary forces immediately begin to affect the flow pattern. 7 The fluid is a 


drawn up to a height equal to the capillary rise above the free wins of the 


Test No. 1 
a © Test No. 2 
No. 3 


1G. —DBASE IEZOMETRIC EAD FOR NCONFINED LOW 
Flow 


water i in the sand. This rise is completed in a relatively short distance. ‘Since 
the | flow surfaces and potential s surfaces must be orthogonal, the potential sur- 


faces will be moved rearward by the u upw ward flow to the capillary yzone. Hence, 
(in which ¥ is the w veight per unit volume and z is the e elevation) v will decrease a 
a and since the elevation z remains constant, the pressure p must decrease . Con- 
sequently, the atmospheric- pressure surface, or what is hereafter r referred t to as 
the free surface, is lowered slightly : as a result of this ‘upward flow. , Once the 


capillary zone is fully developed, the location “ the free surface is not affected 


of the flow, the free surface is a stream suttnes Sale which the fluid within the : 
capillary zone is flowing at less than atmospheric pressure. — In order for the _ 
- fluid to enter the well that is at atmospheric pressure, the pressure in the fluid _ 


be increased it is equal or above pressure by a a decrease 


| 

= 

— 


q 
MULTIPLE WELLS 
in n elevation. ' The resulting flow across the atmospheric-p -pressure re line near the - 
well causes the potential surfaces to lower in order that the potential and stream > 


= 
surfaces remain normal. Hence, at a given point the potential (h = p/y + 2) 


is increased by a an increase in pressure, resulting in a higher atmospheric-pres- 
_ sure line . Thus, the capillary flow has lowered the free surface When the f flow is. 
established at the e periphery, and raised | the free surface in the vicinity y of the 
well, as can be seen from Figs. 4 and 5. 
flow- net +t principles were used to estimate ‘the extent to w hich the free 
surface was altered by the downw ward capillary flow near the Ww well. . Since _— 
_ free surface unaffected by capillary fi flow is a stream surface that must ; be normal — 
“to ‘o each pc potential surface, the atmospheric- -pressure surfaces measured in the 
pote model tests were adjusted until the orthogonality requirement was met, 
as shown i in 1 Fig. 2 The ec correction was considerable and could not be ignored. 
‘When the corrected free surface was is plotted against the innit of the term 


o Test No. a 
@ Test No. 2 
® Test No. 3 


©TestNo.4 | 
6 


0 60 40 
Fra. 7. REE SURFACE as A FuNctTIon oF Rapivs 
(Fig. 7), an essentially straight line was obtait ned, a linear 


function of the 
Gin which m is s the slope and has the dimension ofa length. 
Since a similar equation w was obtained by the writer from a re-examination 
of the Babbitt-Caldwell data, sufficient experimental confirmation of the linear 
my 4 logarithmic nature of the fr ree surface near ‘the well i ‘is available to justify an 
analysis | of the flow pattern to learn the implication of such a relation. 


sider the differential flow (¢ = Q) occurring between the free surface and an 
adjacent. stream surface. _ By assuming that the vertical distance between 
- stream surfaces remains constant (Az = bop C2) as the flow approaches the w vel, 
following linear logarithmic is obtained d by combining the Darey 
and continuity equ: ations, and integrating 


| — 
= 
7 
— 
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Data obtained 1 from a sufficient range of the controlling var are e not avail- 

- 4 able to enable the establishment of the parameters involved in ‘the constants — 
and C2. Additional data should aid materially in defining these t terms. 
However, from Fig. 8, in which m is seni against Q for the single well tests, 7 . : 


can be seen that Q is the 


© Test No. 1 bcs 


primary variable. Hence, the 
assumption the vertical 
“4 ‘distance betw een the free sur- 

face and an adjacent flow sur- 7 
 facer emains constant, not 

results in a form of the equa- 
_ verified experimentally, 


gives a very: good | clue as to 


conclusions regard- 
_ = the nature of the coefficient 
Discharge, Q(x 1 10°) is obtained. When additional 
mination of the “coefficient m 
are available, it is felt that the equation 1 will be more useful if the boundary — 


ia - conditions at the well : are used as the reference rather than those at the radius | 

of influence, as in Eq. 5. 


Under ordinary field is no: definite radius of influence, be- 
cause the conditions far from the well are dependent upon the natural recharge, 
80 an analysis of the variables at the v w vel was undertaken in an attempt vl 
_ define the flow conditions. . WwW hen the variables Q, K, rw, hw, and h, at the w 
are | combined to form dimensionless the f following rela- 


K 


é in which h, is the height. of the intercept of the free surface with the e edge of fhe 
Lag 


celectrical-ar -analogy tests nd the writer’ data from “the anodel tests, the 
parameters w were plotted as s shown in Fig. 9, thereby obtaining the e discharge for 


unconfined well flow in terms of the boundary conditions at the well. — a 


_ of the nature of the dimensionless parameters, a family of curves was quite ac- 
‘curately obtained from the available experimental data. However, such 
presentation fails to give the values of h, , for the ieupertent case of complete 


he being zero for all values of 


— 
| 
oa 
¢ 


MULTIPLE LE WELLS 


ables at the w ataaal were grouped into into the dimensionless is 


“os resulting family ¢ of curve curves (Fig. 10), although less accurate because of the ; 


range of the data, the relationships when when hw - > 


86 — 


=. 


© Babbit- 


Values of of = 


‘Fra. 9. NCONFINED WELL Fiow 1n TERMS oF Bounpary ConDITIONS AT THE WELL 


x To: verify this functional relation further, the values ¢ obtained from the plot 


vells in any arrangement can be defined in terms of the boundary conditions 
at the w ell, and hence, the vague quantity—radius ¢ of influence—can be avoided. 
It ‘should | be remembered that both dimensionless plots are based on limited 

Consequently, the results should be taken as a positive indication of 
possibilities « of st such an analysis, and not as the final answer. WwW hen t the results 


of “other experimental work : are re available, the usefulness of the plots ca can then be 

re ~ Since the free surface next to the well is not only hard t to. measure, but is also 
affected considerably by the entrance conditions, the intercept at the well was 

obtained by plotting the free-s -surface elevations against the logarithm of 


term r/r and extending the essentially - straight line to the well. (Fig. 7). 


similar procedure i is recommended when field data are used to relate the “al 


good agreement. This means that unconfined flow to either s or ‘multiple 


The fundamental nature of the dimensionless parameter Kr should not 
be. overlooked. This p parameter is an index of the shape ¢ of f the cone of depres-_ 


sion ¢ of thet e wa vater table around the Ww ell, large values of the 


| 
| 
> 
: 
— 
& Bi 
| 


>. 


An analysis will show that this shee fun 


= 


damental i in all well flow, 


; is nothing more than the ratio of the Froude and Reynolds n numbers, which nn 


: the ratio of the viscous to the gravity forces; - the inertial forces cancel because ~ - 
their minor importance in | the viscous and forces. 


2 


nN 


Values of ; 


Wey 


10. .—UNCONFINED WELL DISCHARGE AS s A FunNcTION OF DrawDowNn AND ‘Free Sunrace NTE PT 


Kr, 


Gin which is the viscosity) is th 


a, through Ww which the flow occurs, the quotient =—; has ¢ 
ity, that, could be considered as s characteristic | of a 


boundary ec condition. Substituting v’ for 
== ~ OF 
well parameter, ‘the following equation is obtained: 


, together with y = pg, into the 


permeability K= 


e following results” 


ering as a representative area characterising the 


tn 
the dimensions of veloc- 

flow produced by a given 


(146) 


— 
bu e im 
‘ 
| 


MULTIPLE WELLS” 
a in w hich F and R are the Froude and Reynolds numbers, respectively. How- 


- 


_ viscous fore ces” 


to the gravity forces ‘the the cone of 


sion around the well. Since this ratio has as its primary variable the discharge 
of the well, it is san that - be called the discharge number. © 
Effectiveness. common expression in well hydraulics 
: bth, r used as a measure of the losses occurring at the well casing. | 


wan 


7 or +i ratio betwe een the « draw down measured 0 d outside and 1 inside the well. 


- If the piezometric heads in 1 Eq. 17 are measured in the zone in which Dupuit’s Ss 

equation . applies, the results are indicative of the losses occurring at the well — 

casing. However, when the piezometric heads” are measured outside the the 
- Dupuit reg region and ‘in partioular at or near the free e surface, the definition g given | 

is not only erroneous but grossly misleading. Inf fact, from F Fi ig. 9 it is seen t that : 


varies ‘from 0 to 1 0, depending on values of Nwo/Tr a and 


the possible. extremes of ha/he of 0 and 1 1. 0: are re substituted into E q. 17, it is seen 


4 that the value of Ew can be varied over the entire r range ge from 0 to 1.0 without 7 
an Too often these basic | principles are ignored, with the result that the s seepage » 
face that must accompany this type of flow is considaned as well loss caused by x 
a poor design or poor development of the well. ‘ae When s subsequent ext extensive de- _ 7 
: velopment fails to alter the flow pattern materially, the premature aeedinie 
is too often reached either the know his profession or the 


‘trapolated to its the casing. ia losses c occurring gat the 
well casing can then be determined by comparing the extrapolated value of h, 
_ with the value cal calculated | by the use of Fig. 9 9 or Fig. 10. aa The ratio between these 
two values is a measure of well efficiency * When the square of the 5 piezometric : 
4 heads plotted against the logarithm | of the radius neues ina a straight lin line, , the _ 


then 
=F 

— 

| 

oY 

— 
the free surface near the well is observed, the piez j 


— to the de depth of of the Ww ater in) the w w all, , will also be a a measure 0 of w ll. 
_ efficiency. ‘The loss oceurring at the well will be the difference bet een these 


The loss oc occurring at the w ell c casing may b ‘be the result of vn or both of two 
“factors: (1) In the derivation of the well equations, s, the pe permeability w was as 


to be constant throughout the porous | media. silt and clay ¢ are carried 


by the: water to » the casing | and their | presence tends to restrict the flow | passages, 
the permeability will decrease, resulting in an increased loss of energy, and (2) _ 

_ int the derivation, it was assumed that Darcy’s law applied; however , Darcey’ s 
hawa applies « only when the flow is laminar. If the flow becomes turbulent, the 
“Hen head i is proportional to second power of of the 


To clarify further | the fact that the seepage face i is not aw rell loss, the dis- 
charges from unconfined and confined wells of similar. geometry and location - 


“will be compared. hen the. ) Dupuit for unconfined flow is solved 


the confined discharge i is 


The value of represents so some peripheral in so 
& the depth | of water in each well (hw) in addition to K, r, and r, are equal; 


then by dividing the unconfined discharge (Qu) by the confined discharge (Q.) : 


: 
an inter esting equation is is obtained 


of the seepage face does not result i in a less w vell system. 


The potential for unconfined was found, by model experi- 
ments on single and multiple wells, to extend essentially unaltered across the © 
_ atmospheric- “pressure surface and completely through the capillary zone into’ : 
the capillary fringe. In model studies in which the heigl | 
pillary fringe. bs n model studies in which the 1eight of the ‘capillary 1 rise 
is significant fraction of the size, hey of the flow on 


ee. : indication is that the measurements were taken in the Dupuit zone. The a 
a 
“ght 
— 
— 
9 
BY’: 
ee ___ Since ¢ must be less than hy, which must likewise be less than h (¢<hw< a 
the unconfined discharge must always be greater than the confined dis- 
— 
iii 


for a single wi well ‘under various: drawdowns is corrected for - the effect of the 
capillary flow, the shape of the free surface near the well can be closely ap- | 
_proximated by a linear logarithmic ft function. ‘The fi free surface i in unconfined 
well flow must always intersect the well casing above the water surface in ‘the : 
well, resulting i in a zone of seepage that increases as the drawdown i increases. 
‘The free surface, consequently, i is above that calculated by the Dupuit equa- 
tion, the difference being greatest at the edge of the well. papell pen Bis ll 


_ Measurements: the base piezometric head show a good | agreement with 


Based upon experimental data, functional relationships s (Figs. 9 a 


_ between the variables at the well have been developed for unconfined flow. 


These ‘relationships : apply to the solution of either single or multiple systems. 


dir dimensionless parameter characterizes the Shape o of the cone of de- 


pression occurring around a well. This parameter, fundamental in all well 

+e flow, is nothing more than the ratio of the Froude and Reynolds numbers or 


the ratio of the viscous to the gravity forces. Since this ratio has as its 

4 rimary variable the the discharge of the well, it is suggested t that = s~ be ealled 


“the discharge number. 
commonly used term “well effectiveness” is grossly ‘for un- 
 ~iiad flow unless piezometric heads are measured in the zone in which the © 


Dupuit equation applies. © Ww hen the square of the piezometric heads plotted 
: against the logarithm of the radius results in a straight line, the indications a are 
that the measurements were taken in the Dupuit : zone. — The intersection of : 
the resulting straight line, extrapolated to the well c casing, w vhen compared to the 
depth of the water in the well, will be a measure of well efficiency. ‘The loss” 
_ occurring at the well will be the difference between these two values. If the 
4 - surface near r the \ well is observed, the piezometric head should be plotted — 
- against the logarithm ¢ of the radius and the e resulting straight line extrapolated | 
until it intersects the well. casing. ;. The losses 0 occurring near the well casingcan 
then be determined by comparing the extrapolated value of h, with that cal- — 
= from the curves of Fig. 9 or Fig.10. The ratio between these two values” 
a measure of well efficiency. W ell | losses can be attributed to either a 
- decrease in permeability near the well or to turbulent flow. _ The zone of 
“4 seepage i is not a well loss. 


paper is based on a entitled “Evaluation of 


pte curately describe the base piezometric head distribution but also describe the 
_ free surface with increasing accuracy as the radius increases. Over the zone => ; 
/ in which the Dupuit equation applies, the stream-surface spacing is propor- _ — 7 
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a _ = = Babbitt- Caldwell dimensionless, variable coefficient relating the draw- 


. 


down to the radius; 


h piezometric head, p/y + z: 


piezometric head at the radius of influence; 
sau of the = the free s surface with the edge of 


he piezometric head in the w sell; sedans 
coefficient of permeability having the dimension of velocity; 


- is slope of the linear logarithmic equation representing the free surface 


near 


pom radial distance from axis Ww ell: vot 
length along a stream line; wile 
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